Long distance entanglement in one-dimensional quantum systems under sinusoidal 
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We investigate entanglement generation in one-dimensional quantum spin systems with the si- 
nusoidal deformation. In the systems, the energy scale of local Hamiltonian at the position x is 
modified according to the rescaling function sin" [-^ (a; — |)] , where A'^ is the length of the system. 
We show that at zero temperature the system with a > 2 is able to generate a sizable entanglement 
between two spins at open edges even when the two spins are infinitely far apart. This long-distance 
entanglement is rather robust against thermal fiuctuations and survives up to a temperature that 
decays with the system size slowly, in an algebraic form. 
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I. INTRODUCTION 



Entanglement is an essential resource in quantum in- 
formation tasks. [l| In particular, the generation of a large 
entanglement between two parties located far away from 
each other is a crucial ingredient in quantum information 
processing such as quantum computation and quantum 
teleportation. One possible way to create an entangle- 
ment between distant parties is to connect them by a one- 
dimensional (ID) quantum many-body system including 
quantum spin chains and correlated fermion/boson sys- 
tems. Such a quantum system which is able to sustain 
a sizable entanglement between qubits at a distance has 
been sought under the concept of long distance entangle- 
ment (LDE). 

The realization of LDE is a challenging task since en- 
tanglement in a system with short-range interactions usu- 
ally decays quite rapidly. However, the considerable ef- 
forts in recent years have revealed that there are some 
systems which can generate LDE.[2|-[l3| A typical exam- 
ple of the systems exhibiting LDE is the spin- 1/2 bond- 
alternating chain with two additional spins weakly cou- 
pled to the open ends. [3, 3 The ground state of the sys- 
tem consists of entangled pairs of spins (singlet pairs, in 
the words of quantum magnetism) in the bulk and two 
effectively-free spins at the open edges. At a sufficiently 
low temperature, the edge spins form a singlet pair across 
the bulk spins and realize LDE. A similar mechanism ap- 
plies also to other systems with gapful excitations. [1,0-131 
Setups utilizing a ID critical spin system[2l-[l] and bosons 
in a ID optical latticeQ to mediate LDE between edge 
qubits as well as dynamical setups to generate LDEj8l-[l3| 
have also been proposed. 

In this paper, we propose a new class of ID systems 
which exhibit LDE, that is, the systems under the si- 
nusoidal deformation (SD').|14|-[2G| In the systems, the 
energy scale of the local Hamiltonian at a position x is 




FIG. 1: Rescaling function /i"' of the sinusoidal deforma- 
tion, Eq. fl} with a = 1,2, and 3. /i"' is maximum at the 
center of the system, x = (A^ -I- l)/2, and becomes zero at 
open edges, x = 1/2, N -\- 1/2. 



modified according to the rescaling function. 
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where TV is the number of sites in the system. The energy 
scale is a maximum at the center of the system, x = 
(A^-|-l) /2, decreases smoothly as the position x goes away 
from the center, and becomes zero at the open edges, 
X = 1/2, N -\- 1/2 (see Fig. [1]). We show that in the 
ground state of the system with a > 2 the entanglement 
between two spins at the open ends of the system remains 
finite even at the limit of infinitely-long system, that is, 
LDE is realized. In particular, for a = 2 the emergence 
of LDE can be proven exactly. The amount of LDE in 
the ground state increases with a. We also examine how 
robustly the LDE survives at finite temperatures. It is 
found that the LDE becomes more fragile against thermal 
fluctuations as a is larger; In the A^-site system with a > 
2, the temperature at which the entanglement between 
edge spins vanishes decreases with N in the scaling form 
~ N~°'. Our results provide materials for flnding the 
optimal a to maximize the amount of LDE at the lowest 
temperature one can reach. 
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The paper is organized as follows. The definition of 
models considered in this paper is presented in Sec. Ill AI 
In Sec. Ill B| we review known results for the systems 
under SD with a = 2 and show that LDE between edge 
spins is realized in the ground state of the systems. We 
then present our numerical results for the ground state 
and finite temperatures in Sec. IIII AI and IIII Bl respec- 
tively. Section HVl contains the summary and concluding 
remarks. 



II. SINUSOIDAL DEFORMATION 



Model Hamiltonian 



We study the spin-1/2 antiferrornagnetic XXZ chain 
under the sinusoidal deformation. |14H2Cl| The model 
Hamiltonian is given by 



N-l 



1=1 

where Si = {Sf , Sf , Si) is the spin-1/2 operator at Ith 
site and A is the parameter of the exchange anisotropy. 
The rescaling function /i"-* is defined in Eq. ([1]). We 
set the exchange coupling constant J = 1 throughout 
this paper and consider the cases of the Heisenberg chain 
(A = 1) and the XY chain (A = 0). Note that the spin- 
1/2 XY chain is mapped to the ID system of free spinless 
fermions by Jordan- Wigner transformation. 



B. Long-Distance Entanglement by Sine-Square 
Deformation 



The SD was originally introduced in Ref. Il4l as an im- 
proved version of smooth boundary conditionj2ll [2^ to 
suppress open boundary effects. In particular, for the SD 
with a = 2, which is called the sine-square deformation, 
an interesting phenomenon was found: Despite of the 
presence of the open boundaries, the boundary oscilla- 
tions in local quantities were suppressed completely. Q 
Subsequently, it was found that a ID critical (gapless) 
system with the sine-square deformation had the ground 
state which was identical in the level of the wave func- 
tion to the one of the same uniform system under the 
periodic boundary condition, [l^ [2^ The equivalence of 
the ground-state wave functions was proven rigorously for 
several systems including ID free-fcrmion system [equiva- 
lent to the XY chain, Eq. ([2]) with A = 0] , the transverse- 
field Ising model at criticality, and the Gaussian model of 
the c = 1 conformal-field theorv ^l6^,17. . .20| . Furthermore, 
it was found numerically that the same phenomenon oc- 
curred in a wide variety of ID models with gapless ex- 
citations such as the spin-1/2 XXZ chain and two- leg 
ladder in magnetic field 11 5| and the (extended) Hubbard 
chain in a metallic phase [18|. These results suggest that 



the equivalence of the ground state between systems un- 
der the periodic boundary condition and the sine-square 
deformation is a generic feature of ID critical systems be- 
longing to the universality class of Tomonaga-Luttinger 
liquid. 

In relation to LDE, the above observation immediately 
leads to an intriguing conclusion: In the ground state of 
the system with the sine-square deformation, the spins 
at the open edges, which are at a distance of the length 
of the system, behave in completely the same way as 
the nearest-neighboring spins in a periodic system. Since 
the entanglement between the neighboring spins in the 
periodic chain remains finite even in the thermodynamic 
limit, it follows that the edge spins in an open chain with 
the sine-square deformation also sustain the same finite 
amount of entanglement. It is thereby concluded that 
the system under the sine-square deformation exhibits 
the LDE. We wiU show later that the model Q for a > 2 
also realizes the LDE in the ground state. 



III. NUMERICAL RESULTS 

In this section, we present our numerical results of the 
entanglement between edge spins in the spin chain under 
the sinusoidal deformation. To evaluate the entangle- 
ment between two end spins, we calculate the concur- 
rence defined as follows. [23. [25j 

We start from the density matrix (DM) for the whole 
system, i.e., p = |$o)('I'o| at T = (|$o) is the ground 
state) and p = for finite temperatures T > (Z is 

the partition function). We set the Boltzmann constant 
kB = 1- The reduced DM for the edge spins Si and Sat is 
obtained as piN = Trj^p 
for the bulk spins at Z = 2, — 1 
for the edge spins is then defined as 



where Trjjy denotes the trace 
The concurrence 



Cip 



IN) 



max(0, Ai — A2 — A3 — A4) 



(3) 



where Ai (Ai > A2 > A3 > A4) are the square roots of the 
eigenvalues of the non-Hermitian matrix PinPin- Here, 
Pin, called the spin-fiipped state, is defined by 



Pin = {cr\ ® CT](f)pi^(crj' (g) cr^). 



(4) 



where cr^ is the y component of the Pauli matrix. The 
concurrence is equal to 1 for the maximally entangled 
state while it is for a separable state. 

Using relations between the matrix elements of the re- 
duced DM Pin and expectation values of edge-spin op- 
erators, a useful expression of the concurrence has been 
obtained, [2^ 
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FIG. 2: End-to-end concurrence C(pijv) in the ground state 
of the system with the sinusoidal deformation [Eq. ((2|] as 
functions of system size N; (a) Heisenberg case (A — 1) and 
(b) XY case (A = 0). Estimates of C{pin) ai, N ^ oo are 
shown in (c) as functions of a. Lines are guides for the eye. 



concurrence for the Heisenberg (A = 1) and XY (A = 0) 
cases, which exhibit essentially the same behavior. It is 
clear that for a > 2 the concurrence takes a finite value 
even at the thermodynamic limit N oo. For a = 4, we 
could not obtain the data for > 28 since the DMRG 
calculation did not converge the true ground state due 
to the extremely small energy scale around open edges. 
However, the exact-diagonalization results for < 24 are 
enough to suggest the convergence of the data to a finite 
large value at iV — >■ oo. For a < 2, on the other hand, 
the end-to-end concurrence decreases slowly, presumably 
in a power law, with the system size N. We note that 
the concurrence for a = 1.9 decreases very slowly but 
monotonically with N. From the results, we conclude 
that the system with a >2 realizes the LDE. 

Figure EJ^c) shows the end-to-end concurrence in the 
thermodynamic limit for the LDE regime as a function 
of a. As estimates of the concurrence at A^ — > oo, 
C(pijv) for N = 64 (N ^ 24) is plotted for a < 3.0 
(a = 4.0). The figure suggests that the end-to-end con- 
currence in the ground state is monotonically increasing 
with a. This is consistent with an intuitive argument 
that the correlation between an edge spin Si (S^r) and 
the neighboring spin S2 (Sn-i) becomes weaker as the 
ratio of the edge bond to the bond next to the edge is 
smaller, which results in a larger entanglement between 
the edge spins from the monogamy condition. We note 
that for a = 2 the end-to-end concurrence at T = 
is equivalent to the concurrence between the nearest- 
neighboring spins in a uniform periodic system as dis- 
cussed in Sec. Ill Bl Indeed, the numerical results for 
a = 2 coincide with the exact value obtained from the 
nearest-neighbor spin correlations in the infinite uniform 
spin chain: limAr^oo C(piAr) = 0.3863... for the Heisen- 
berg chain (A = 1) and limAr_>oo C(piAr) = 0.3393... for 
the XY chain (A = 0). 



B. Finite temperatures 



where M = (Sf) = (Sf^) and Cfjv = {SfS'^) {a = x,z). 
We calculate numerically the correlation functions CJjy 
and C^j^ at zero and finite temperatures (note that M ~ 
in our model without external magnetic field) , and then 
evaluate the concurrence using Eq. ([5]). 

A. Ground state 

In Fig. [21 we show our numerical data of the end-to- 
end concurrence in the ground state (T = 0) of the SD 
systems with up to A^ = 64 spins. The data for A^ < 24 
were obtained by the exact-diagonalization method while 
those for N > 28 were calculated by the de nsity -matrix 
renormalization group (DMRG) method. [13, Figures 
[5] (a) and (b) show the system-size dependence of the 



In the previous section, we have shown that the SD 
system with a > 2 exhibits the LDE. The end-to-end 
concurrence in the ground state increases with a, sug- 
gesting that larger a is preferable for achieving a larger 
amount of LDE. However, since the SD with large a leads 
to a small energy scale around the open edges, it is nat- 
urally expected that the system with large a becomes 
fragile against thermal fluctuations. Therefore, it is im- 
portant to examine how robustly the LDE in the system 
with SD survives at finite temperatures. To the end, we 
have calculated the end-to-end concurrence at finite tem- 
peratures by the quantum Monte-Carlo method based on 
the directed-loop (worm) algorithm. [29| - [3T| 

In Fig. [31 we present the numerical data of the end-to- 
end concurrence at finite temperatures for A^ = 24 as a 
typical example. For all the cases of a and A calculated, 
we have observed qualitatively the same behavior: When 
T increases from zero, the concurrence remains at almost 
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FIG. 3: Temperature dependence of the end-to-end concur- 
rence C(pijv) in the system under the sinusoidal deformation 
[Eq. ((2|] with A'^ = 24 spins; (a) Heisenberg case (A — 1) and 
(b) XY case (A = 0). All data are plotted with error bars 
but the most of them are smaller than the symbol size. Lines 
are guides for the eye. 



the same value as that for T = for some temperature 
range, and then decreases smoothly down to zero. As a 
is larger, the concurrence is larger at T — > but starts to 
decrease at a lower temperature, as expected. [12] On the 
other hand, when a is small, the concurrence is rather 
small at T but can survive up to high temperatures. 
Therefore, there is a trade-off between the amount of 
LDE achievable at T = and the temperature range for 
which the LDE is sustainable. For a practical purpose, 
one must choose a with considering how large LDE is 
desired as well as how low temperatures one can reach. 
Our data may serve as a basis for the judgement. 

To evaluate the robustness of LDE more quantitatively, 
we examine the temperature T* at which the end-to-end 
concurrence C(piAr) vanishes. Figures|3](a) and (b) show 
the system-size dependence of T* . The data clearly sug- 
gest that T* decreases with TV algebraically, T* ~ . 
The decay exponent 77 obtained from the fitting is plot- 
ted in Fig. 131(c) as a function of a. For a > 2, the expo- 
nent obeys a relation = a. This behavior is not triv- 
ial but can be understood from the fact that the energy 
scale at edge bonds, which should give the energy scale 
of the excited state localized around the edges, changes 
as sin"(7r/A'') ~ 7V~". For a < 2, on the other hand, the 
decay exponent 77 is larger than a. While this numerical 
observation indicates the presence of low-energy excita- 
tions which has an excitation energy smaller than that of 



10 20 30 40 50 60 70 




2 3 
(C) « 



FIG. 4: Temperature T* at which the end-to-end concur- 
rence vanishes as functions of system size A'^; (a) Heisenberg 
case (A = 1) and (b) XY case (A = 0). Lines show the fits 
to an algebraically decaying form T* = ^A'^"'' with fitting 
parameters r; and A. (c) Decay exponent 77 as functions of a. 
The dotted line represents the relation rj = a. All data are 
plotted with error bars. 



the edge states, it is not clear how to construct such an 
excited state. This question for a < 2 is open for future 
research. 



IV. CONCLUDING REMARKS 

We have studied the long distance entanglement be- 
tween edge spins of onc-dimensional spin systems under 
the SD. fn the systems, the energy scale of the local 
Hamiltonian is modified smoothly from a maximum at 
the center of the system to zero at open edges according 
to a rescaling function Eq. ([1]). When the exponent in the 
rcscaling function is a = 2, it was shown that the ground 
state of the system was equivalent to the one of the cor- 
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responding uniform periodic system. Therefore, the sys- 
tem under SD with a = 2 generates an entanglement 
between edge spins whose strength is the same as that 
between nearest-neighboring spins in the periodic chain. 
The system thereby reahzes the long-distance entangle- 
ment, a finite entanglement between edge spins infinitely 
far apart from each other, in the ground state. 

We have investigated numerically the entanglement be- 
tween the edge spins in the SD system for various a. Us- 
ing the exact diagonalization and DMRG methods, we 
have calculated the end-to-end concurrence in the ground 
state. It is then found that the SD with a > 2 gen- 
erates LDE. The amount of the entanglement is larger 
as a is larger. We have also examined how robust the 
LDE is against thermal fluctuations. The numerical data 
of the cnd-to-cnd concurrence at finite temperatures ob- 
tained from the quantum Monte-Carlo method suggest 
that when temperature increases, the concurrence keeps 
the value at T = for some range of low temperatures, 
then disappears smoothly. The temperature T* at which 
LDE vanishes decays with the system size TV in a power 
law, T* ^ N~^. For a > 2, the decay exponent rj is 
found to be equal to a. Our numerical results thus in- 
dicate that large a is preferable for generating a large 
LDE at T = but disadvantageous in sustaining LDE at 
finite temperatures. For a practical purpose, one may se- 
lect the optimal a in the energy deformation considering 
the amount of LDE required and the lowest temperature 
which can be achieved. 

A characteristic feature of SD is that it realizes a "true" 
LDE in the ground state in the sense that the entangle- 
ment between edge spins remains finite even at the in- 



finitely long chain. Furthermore, the LDE can survive 
up to a rather high temperature; the temperature T* at 
which LDE vanishes decays quite slowly, in a power law, 
with the system size. The latter point is in contrast with 
the LDE mediated by systems with gapful excitations, for 
which T* typically decays exponentially with system size. 
We also emphasize that our results are not restricted to 
the spin-1/2 chains but applicable to general ID critical 
systems belonging to the universality class of Tomonaga- 
Luttinger liquid. These properties of SD would make it 
be of wide application in entanglement engineering using 
quantum correlated systems. 

The realization of a system with SD in laboratory is 
an intriguing and challenging problem. The systems of 
ultra-cold atoms in optical lattices, for which hopping 
amplitudes and interatomic couplings can be tuned [33|- 
[35| . should be a promising candidate. 
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